Scale- free antiferromagnetic fluctuations in the 5 = 1/2 kagome antiferromagnet 

herbertsmithite. 
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Neutron spectroscopy and diffuse neutron scattering on herbertsmithite [ZnCu3(OH)6Cl2], a near- 
ideal reahsation of the s = 1/2 kagome antiferromagnet, reveal the hallmark property of a quantum 
spin liquid; instantaneous short-ranged antiferromagnetic correlations in the absence of a time- 
averaged ordered moment. These dynamic antiferromagnetic correlations are weakly dependent of 
neutron-energy transfer and temperature, and persist up to 25 meV and 120 K. At low energy 
transfers a shift of the magnetic scattering to low Q is observed with increasing temperature, pro- 
viding evidence of gap less spinous. It is argued that these observations provide important evidence 
in favour of RVB theories of (doped) Mott insulators. 
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There has been a long search for quantum spin liq- 
uids in Mott insulators in which the common Neel anti- 
ferromagnetic order is destabilised by quantum fluctua- 
tions [1, 2'], the mixing-in of spin-singlets. The 5 = 1/2 
kagome antiferromagnet is one system where such a spin 
liquid state could, in theory, be found. The kagome lat- 
tice, which owes its name to a Japanese basket weaving 
method, is a net of corner-sharing triangles. When anti- 
ferromagnetically (AF)-coupled spins are arranged at the 
vertices of this lattice it becomes impossible to satisfy all 
AF bonds simultaneously. Due to this geometric frustra- 
tion, the ground state has a continuous and macroscopic 
degeneracy. 

Both in the classical limit of large spin [5, and in 
the quantum limit of s — 1/2 [5], the kagome AF has 
been predicted to retain the full symmetry of the mag- 
netic Hamiltonian. Exact diagonalization (ED) results 
suggest that in the macroscopic limit of the 5 = 1/2 
system a continuum of non-magnetic global spin-singlet 
<5tot = states fills an energy gap to magnetic excita- 
tions . Whether there is a spin gap in real systems has 
however been debated and in most theories building on 
Anderson's RVB proposal [2j], such as the algebraic spin 
liquid (ASL) [7, 8, 9], there are gapless spinous [tI. fiot. 
Spinous are 5 = 1/2 excitations which are created in 
pairs following a singlet-triplet excitation but can disso- 
ciate with no (additional) cost in energy. The ASL is 
thought of as a "mother of competing orders" [11] in- 
cluding Neel antiferromagnetism and d-wave supercon- 
ductivity. As the insulating analogue for the Fermi liq- 
uid in metals the ASL would be a natural candidate for 
the ground state in a Mott insulator where Neel order 
is suppressed. The results presented here on the kagome 
AF compound herbertsmithite provide further evidence 
for this idea. At the same time, the absence of a spin gap 



in herbertsmithite can be reconciled with ED results [12| 
when a Dzyaloshinky-Moriya interaction [13] or struc- 
tural disorder is taken into account. 



Of the many kagome systems 
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mentally, including the jarosites 
(SrCr8-a:Ga4+a;Oi9) and volborthite 
smithite [ZnCu3(OH)6Cl2] is the first where no (par- 
tial) freezing of the spins is observed down to the low- 
est temperatures accessible experimentally [1^,0. The 
Weiss temperature is 300 K and the AF exchange inter- 
action J ^V^^Yi [6*, ^20*]. Herbertsmithite contains well 
separated 2D kagome layers of Cu^+ ions linked by 
ions in 0H~ groups. Separating the kagome layers are 
Zn sites of Oh symmetry, which can also host Cu^+ ions 
to form the zinc paratacamite family ZnajCu4_aj(OH)6Cl2 
with < X < 1. For x = 1 the low-temperature suscep- 
tibility is dominated by a Curie-like contribution from 
"antisite spins" (a ~ 6% fraction of Cu^+ spins which 
have traded places with Zn^+ to occupy the interplane 
Zn site) [2l|, 



22[. These weakly-coupled 5 = 1/2 spins, 
which have been identified as individual doublets [2l|, 
have also been observed in neutron spectra at the Zee- 
man energy in applied fields [13, [l^. Using ^^O NMR [24 ] 
and CI NMR |25| there is now convincing evidence that 
the kagome layers have a nonzero susceptibility as T ^ 0. 
This is in agreement with the continuum of magnetic ex- 
citations observed using neutron spectroscopy at energies 
between ~ 0.8 and 2 meV [2C|. Moreover, this spectrum 
is independent of the temperature, which has been sug- 
gested to point to the proximity of a quantum critical 
point [20]. Here we present neutron spectroscopy data 
up to energy transfers of ~ 30 meV and temperatures 
up to 120 K, to measure the dynamic magnetic structure 
factor and fully characterize the quantum spin dynamics 
in herbertsmithite. 
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Large single crystals of herbertsmithite do not 
presently exist. Hence, 20 g of 98% deuterated herbert- 
smithite powder was synthesised following the method 
as described in Ref. 18. This sample was character- 
ized using neutron diffraction, DC magnetic susceptibil- 
ity, heat capacity [2l| and /iSR measurements [19]. The 
measurements were carried out using the polarised neu- 
tron spectrometers D7 and IN22, and the time-of- flight 
(TOF) spectrometers IN4, ILL, France, and MARI, ISIS, 
UK. At D7 the magnetic scattering was separated from 
the nuclear and spin-incoherent scattered neutrons using 
XY Z polarisation analysis [2^. No energy analysis was 
carried out on the scattered neutrons, with an incident 
neutron energy of 8.95 meV, effectively energy integrat- 
ing the neutron cross section [27] up to energy transfers 
(e) of ~ 6.5 meV. Complete S{Q^e) maps without po- 
larisation analysis were obtained at 2, 4, 10, 30, 60 and 
120 K at IN4 using 17.21 meV neutrons and at 2K at 
MARI using 56 meV neutrons. The normalised spin and 
nuclear incoherent signals measured at D7 were used to 
normalise the IN4 and MARI data. Neutron polarisation 
and energy resolved measurements were carried out at 
the triple- axis spectrometer IN22, ILL. Several Q, e and 
T scans were performed for spin flip (SF) and non spin 
flip (NSF) channels for incident and scattered neutron 
spins polarised parallel to Q (labelled xx) and perpen- 
dicular to Q (labelled zz). 
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FIG. 1: (Color online) a) Instantaneous magnetic correlations 
at 4 K (open circles) and 10 K (filled circles) from D7. Data 
taken at 60 K (not shown) resembles the 10 K data. The 
solid lines are a guide to the eye. b) The Q dependence in 
the dynamic correlations from IN4 and MARI with the en- 
ergy integration interval indicated in the legend. The dotted 
line in panel a and b is the structure factor for dimer-like AF 
correlations, for the dashed line a single-ion contribution cor- 
responding to the 6% antisite spins in this system is added, 
c) The energy and temperature dependence at Q = 1.3 
from the TOF data. 

At D7 the energy of the scattered neutrons is not anal- 
ysed. The D7 data (Fig. [U^a)) therefore correspond to 
magnetic correlations which persist on a time scale of at 



most [23 - 6.5 meV (1.6 THz). The broad peak ob- 
served here at 1.3 is direct evidence of short-ranged 
instantaneous near- neighbour AF correlations. For com- 
parison the powder- averaged structure factor for uncor- 
related near-neighbour AF dimers given by 
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is shown in Fig.[T](a) and (b). In Eq.[T]F((5) is the Cu^+ 
magnetic form factor and d = a/2 = 3.42 A the Cu-Cu 
distance within the kagome plane. 
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FIG. 2: (Color online) Inelastic intensity for 3 polarisation 
channels: NSF^x contains phonons and NSF background. 
SFxx contains 2/3 of the magnetic intensity and SF back- 
ground. SF^;^ contains 1/3 of the magnetic intensity. Hence, 
6'mag/3 = SFxx— SFzz- It is seen that magnetic intensity at 
the maximum at Q = 1.3 (SF^x in panel a) is finite and 
nearly constant for all covered temperatures (b) and energies 
(c,d). Data in (a),(c) and (d) were measured at 2K. 

Similar short-ranged AF correlations are observed in 
the TOF data from IN4 and MARI at all e accessed, up to 
25 meV. Fig. WO^) and (c) show respectively the momen- 
tum and energy dependence of the magnetic scattering 
peak in the TOF data. Fig. [2] shows scans in e, temper- 
ature and Q using polarisation analysis which confirm 
that these dynamic correlations are of magnetic origin. 
That a similar Q dependence is observed in the D7 data 
suggests that the dynamic magnetic correlations remain 
unchanged down to the elastic line. From Fig.[2fb,c,d) it 
is clear that in addition to magnetic inelastic scattering 
('S'mag) there is significant inelastic background from spin 
incoherent (SFinc = SF^^ — ^mag) and nuclear scattering. 
This is ascribed to incoherent scattering from residual 
protons and multiple scattering. This background is also 
present in the IN4 data and to a lesser extent in the 
MARI data because there a smaller sample was mounted 
in an annular geometry. The integrated intensity of the 
dynamic AF correlations as found in Fig. [1] and [2] up to 
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25 meV corresponds to 30 to 50% of the sum rule 5(5 + 1) 
for the spins on the kagome lattice. 

In the raw IN4 data the dynamic correlations at Q < 
1.8 A-i persist up to 30 K while the magnetic scattering 
at the maximum of the peak at 1.3 changes very 
little up to 120 K (Fig. □(c) and Fig.[2](b)). The inten- 
sity at the energy loss side (e < 0) in the TOF data 
obeys detailed balance at all temperatures. This com- 
bined with the magnetic scattering intensity and tem- 
perature dependence in the D7 data implies that at the 
elastic line there must be increased magnetic scattering, 
but still with the same Q dependence. These (quasi) 
static magnetic correlations are reduced with a shift of 
intensity to lower Q as the temperature is increased, in 
what looks like spinon excitations in the D7 data at 10 K 
(Fig. ma)). 

The above observations complement previous neutron 
spectroscopy results on herbertsmithite [lO] , which show 
the inelastic magnetic scattering cross section is energy 
independent between 0.8 and 2 meV apart from a weak 
field-dependent peak which is due to the weakly-coupled 
Cu spins on Zn sites. The Q dependence in this energy 
range is not peaked at 1.3 which is due to the single 
ion contribution from antisite spins. This contribution 
also needs to be added to Eq. [T] to fit the D7 data. The 
absence of a peak around the edge of the first Brillouin 
zone in Ref. |20| down to 35 mK implies there is no signif- 
icant increase of the dynamic magnetic correlations over 
three orders of magnitude in temperature! 

At Q > 1.8 and for T > 30 K the inelastic scat- 
tering is increasingly dominated by phonons. The tem- 
perature dependence of phonon scattering can in general 
be calculated from a temperature independent dynamic 
susceptibility using the fiuctuation dissipation theorem. 
The temperature independence of the magnetic scatter- 
ing observed at low Q prompted us to fit the tempera- 
ture dependence of each pixel in S{Q^ e; T) (Fig. [31(a) and 
(b) for T = 2 K and 120 K respectively) as the sum of 
a temperature independent component Sxy{Q-,^) and a 
component following linear response x"{Qi^)i 
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The 6'AF(Q,e) and x"{Qi^) resulting from the fit are 
shown in Fig. [31(c) and (d) respectively. As expected, 
Sa¥{Q-,^) corresponds to the raw data from IN4 at low 
Q and T < 60 K and there is a good overall agree- 
ment with the structure factor of Eq.[T](Fig.[4l(a)) added 
to a constant background. Phonon dispersion originat- 
ing from the nuclear Bragg peaks is on the other hand 
clearly visible in x"{Q->^)' The excellent fit for all data 
at 6 different temperatures, as illustrated in Fig. [41(b) 
for a small sample of points, confirms that the mag- 
netic correlation length does not diverge as the tem- 
perature is lowered and that these dynamic correlations 
persist up to at least 120 K. The maximum in x" and 
S{Q^ e) around 7 meV that extends to Q < 1.8 A~^ gives 
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FIG. 3: (Color online) The scattering cross section S{Q,e) of 
at 2 K (a) and 120 K (b). The intensity scale is mBarn st.~^ 
meV~^ form. unit~^. The elastic powder diffraction pattern 
is also shown in (b). Panel (c) and (d) show respectively 5'af 
and x'' • The white band at low Q is a gap between detectors. 
Some increase in intensity in 5'af is observed at higher ener- 
gies but this is likely due to the direct beam and larger error 
bars. As shown in Fig. [IJa) the increase in intensity is only 
very small. 



the impression of an increase in magnetic scattering at 
(Q, e) = (1.3 A~^, 7 meV). This signal in at low Q is 
most likely due to some weak multiple scattering at IN4 
because it was not observed in the MARI and IN22 data. 
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FIG. 4: (Color online) a) The Q dependence in 5'af(Q, e) for a 
number of energy transfers, b) The temperature dependence 
for a number of points in S{Q,e = 4 meV) and fits to the 
data using Eq. [2 



Hence, we found that 1) At low frequencies (e < 
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2 meV) there is a shift of intensity to lower Q as the 
temperature increases, which could be due to (gapless) 
spinous. Spinous then also account for the non-zero mag- 
netic susceptibility as measured with ^^O NMR [24]. 2) 
There are dynamic AF correlations with a very weak tem- 
perature and energy dependence, extending over more 
than three orders of magnitude in temperature and two 
orders of magnitude in energy. Hence, the dynamic AF 
correlations are approximately scale free. There is no 
characteristic energy scale in the system, not even one 
determined by temperature as in systems with uo/T scal- 
ing (with u = e/h) and as expected for the ASL 
Although with neutron scattering any long range or- 
der of spin-singlet dimers cannot be detected directly, 
the absence of a spin gap and the temperature indepen- 
dent dynamic correlations observed here and in [20] im- 
ply the absence of a static dimerisation, at least down 
to 35 mK. Short ranged dynamic correlations are com- 
monly observed in quantum spin liquids at temperatures 
well above their respective ordering and spin-glass tran- 
sitions The temperature independence of the 
dynamic AF correlations observed here suggests that, in- 
triguingly, the short-ranged dynamic correlations could 
be a property of the ground state in herbertsmithite. 

A system of particular relevance in the present con- 
text is the spin-glass phase in La2-^Sr.7;Cu04 (LSCO) 
with doping around x = 0.05 [30, ISjJ. At this level 
of doping, just enough for the Neel order to be com- 
pletely suppressed, the correlation length of the dynamic 
AF correlations becomes temperature independent not 
far below the Neel temperature of the parent compound. 
Such a strong likeness with the dynamic magnetic cor- 
relations in herbertsmithite is in agreement with RVB 
theories, which predict that the Mott insulator in the ab- 
sence of Neel order and the pseudogap phases are closely 
related states. One proposal is that both are described 
by a staggered-flux, or algebraic spin liquid , a state 
with algebraically diverging correlations. Here the agree- 
ment with the spatial and temporal correlations in her- 
bertsmithite, and with the spatial correlations in LSCO, 
breaks down. This could be due to weak structural dis- 
order in the respective Cu-0 planes which both systems 
are known to suffer from. To test this possibility new 
herbertsmithite samples are needed with reduced antisite 
disorder. The energy spectrum and powder- averaged Q 
dependence measured here are, apart from the absence of 
a spin gap, also in rough agreement with ED results [12] 
but large single crystals and theoretical predictions on 
the temperature dependence would be needed for a more 
detailed comparison. 
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